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Abstract: Calculation of amplitudes in perturbative quantum field theory involve 
large loop integrals. The complexity of those integrals, in combination with the 
large number of Feynman diagrams, make the calculations very difficult. Reduction 
methods proved to be very helpful, lowering the number of integrals that need to 
be actually calculated. Especially, the reduction at the integrand level technique, 
improves the speed and set-up of these calculations. In this article we demonstrate, 
by counting the numbers of tensor structures and independent coefficients, how to 
write such relations at the integrand level for one— and two— loop amplitudes. We 
clarify their connection to the so-called spurious terms at one loop and discuss their 
structure in the two— loop case. This method is also applicable to higher loops, and 
the results obtained apply to both planar and non-planar diagrams. 
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1. Introduction 

Modern colliders such as the Large Hadron Collider (LHC) (and the Tevatron be- 
fore it was shut down) produce a large amount of experimental data. In order 
to understand the output of these experiments, comparison between very precise 
theoretical results and experimental results is needed. It is clear, from the theoreti- 
cal point of view, that Next-to-Leading-Order (NLO) and Next-to-Next-to-Leading- 
Order (NNLO) calculations with many external legs have to be considered [1]. This 
implies that (very) large loop integrals have to be computed for very many Feynman 
diagrams, which has widely been considered the bottleneck of such calculations. 

Reduction techniques form a way out. The idea of reducing Feynman integrals 
with a large number of denominators to a set of simpler integrals (i.e. with fewer 
denominators) at one loop goes surprisingly many years back [2, 3]. A typical integral 
with n such denominators, in d space-time dimensions, is given by 

f d d q \ 

J q D 1 D 2 ...D n 
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where Di = {q + Pi) 2 — ttl 2 is the denominator of a generic propagator. In [2] the 
authors reduce a triangle (integral with 3 denominators) to bubbles (2 denominators) 
in 2 dimensions while in [3] a pentagon (5 denominators) is reduced to boxes (4 
denominators) in 4 dimensions. We see that the result of the reduction depends on 
the dimension. However, the methods we will use can be applied to all dimensions. 
Our main interest is of course the case d — 4. 

Since a few decades now [4], it is a very well known fact that a generic one- loop 
amplitude is decomposable in terms of scalar integrals, with one, two, three and 
four external legs (in d — 4). Passarino and Veltman [5] used Lorentz invariance 
to express tensor one-loop n— point integrals in terms of m— point scalar integrals 
(m < n). As a consequence, only the evaluation of scalar integrals (integrals with 
trivial numerators) is needed in order to perform a one-loop calculation. 

In another attempt [6] a pentagon-to-boxes decomposition is performed in 4 di- 
mensions. The importance of this paper is that it provides a basis in four-dimensional 
momentum space (the so called van Neerven-Vermaseren basis), which proved useful 
for understanding one-loop reduction. Another important fact about this paper is 
the use of what we call nowadays spurious terms to decompose a scalar pentagon 
to boxes. Spurious terms are terms that, by construction, vanish upon integration. 
Their role will be explained later when we consider reductions at the integrand level 
and we shall see why one cannot avoid them. 

The next big step comes from unitarity methods [7-11]. Instead of working 
with specific Feynman diagrams these methods have a big advantage in that they 
try to decompose the whole one-loop amplitude in terms of the scalar integrals. By 
cutting propagators 1 the rational coefficients of loop integrals are given in terms of 
products of tree amplitudes. In generalized unitarity methods [12-17], the notion 
of multiple cuts is introduced. One can cut more than one propagator to find these 
coefficients. Note that, for d = 4, cutting four propagators essentially determines the 
loop momentum (there is, in general, more than a single solution since the D's are 
quadratic in the loop momentum). 

The Ossola-Papadopoulos-Pittau (OPP) method [18-22] comes as a natural com- 
bination of all the above. Since every integral can be decomposed to scalar integrals 
with up to four denominators (for d — 4), every one- loop amplitude is written in 
terms of coefficients that multiply these scalar integrals. The OPP method works at 
the integrand level [23, 24], which means that for these decompositions to be possible 
one must also include spurious terms. Then one has to find a way to calculate the 
coefficients of the reduction and multiply them with the appropriate scalar integrals, 
using one of the packages available for the evaluation of them(i.e [25, 26]). Finding 
the coefficients is a purely algebraic problem. The method is suitable for a fully 

1 'Cutting' a propagator means that loop momenta are chosen for which one or more of the D's 
vanish so that the integrand becomes singular. One also speaks of 'putting propagators on-shell'. 
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numerical implementation. The OPP method has been widely used so far in many 
one loop calculations (see for example [27-39]). 

As we noticed before, in the case of one-loop calculations a basis for any inte- 
gral is known in advance. Any one-loop integral can be written in terms of scalar 
boxes, triangles, bubbles and tadpoles. However, in the case of higher-loop integrals 
the situation is different. A basis is not known a priori. It is believed that uni- 
tarity methods can also be applied in that case and there are some recent papers 
in that direction, performing decomposition a la OPP [40-43], or using generalized 
unitarity [44-46] at two loops. 

Two remarks are in order here. The first is that the basis of two-loop integrals 
does not include only scalar integrals. It includes integrals that also have irreducible 
scalar products (ISP) as numerators (to some power) that cannot be rewritten as 
existing denominators of the integral. In the one-loop case these ISP are always 
spurious and integrate to zero, but for higher loops this does no longer hold. The 
second remark is that if one is interested in constructing a unitarity-like basis, the 
set of integrals that ends up with is not necessarily a minimal one: the integrals are 
not by default Master Integrals (MI). There might be smaller sets of true MI and at 
two or more loops one can find them by using integration-by-parts (IBP) identities 
[47-49]. 

In this article we prove that any two-loop integral can be written as a linear 
combination of integrals with at most 2d denominators. From this set of integrals, 
one can in principle end up with true master integrals with numerators containing 
ISP 's. The layout of this paper is the following. We start with definitions and 
reduction at one loop. We see why one can have unitarity-like bases at the integrand 
level by writing first the numerator of our integrals (for scalar cases the number one) 
as a sum of coefficients times denominators. We investigate, using simple algebra, 
what is the form of these coefficients in order for our polynomial equation to have 
solutions. Then we repeat the same procedure in the case of two-loop integrands. 
Our method is not exactly the OPP one, but the connection of the two methods will 
become clear. 

2. Reduction at one loop 

2.1 Introduction : reduction with trivial coefficients 

We start with integrands of any given one-loop amplitude. These integrands consist 
of the sum of integrands coming from the Feynman diagrams that contribute to a 
given process and share the same topology; the advantage is that we perform the 
decomposition once instead of reducing every single diagram separately. For that 
reason we deal with integrand-graphs, or iGraphs instead of Feynman diagrams. We 
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give an example of an iGraph of order 5 (pentagon) below, where j = 1 ... 5, 

Dj = D(q + pj) = (q + pj) 2 - m 2 = q 2 + 2(pj ■ q) + ^ , fij = p 2 - m 2 ,(2.1) 




P3 



A one-loop pentagon iGraph. The dots serve to distinguish the 
several denominators. The external momenta are indicated. 

The loop momentum is denoted by g M , and p^ is called the external momentum, 
where it must be realized that by this we do not mean a momentum related to a 
particle incoming or outgoing in a given amplitude: what we call external momenta 
are simply fixed momenta, given in some way by the configuration of incoming and 
outgoing momenta and the various diagram topologies. 

Consider a one- loop iGraph of order n: 

1 

D 1 D 2 D 3 D 4 ...D n • 

We say that we can decompose this iGraph if we can find funtions Ti j2j ... jn (g) such 
that 

D 1 T 1 (q) + D 2 T 2 (q) + ■■■ + D n T n (q) = 1 , (2.2) 

for then we have 

1 = Uq) T 2 (q) T n (q) 

D 1 ---D n D 2 D 3 D A ...D n D x D 3 D A ...D n D l D 2 D 3 ...D n _ l ' l ' J 

and the original iGraph is decomposed into a sum of iGraphs of order n — 1 (or 
lower) . 

This immediately leads us to state the following theorem: one-loop iGraphs of 
order d or smaller cannot be decomposed in the above manner. The reason is simple: 
for n < d there exist a cut through all propagators, so that Dj = for j = 1, . . . , n 
and eq. ( |2.2| ) then would become = 1. 

The simplest possibility for the functions Tj(q) is to take them to be just numbers 
independent of g M ('trivial' coefficients): 

Tj(q) = Xj . 
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From eq. (|2.2|) then we have 



<r J2 x o + '-'//' Yl •'•'/''" + Yl x sN = 1 • ( 2 - 4 ) 

j=l j=l j=l 

Since this has to hold for any value of g M we must have separately 

n n 

5> i = , = , (2.5) 

3=1 3=1 

and 

n 

X] •'•,/', 1 • (2-6) 

3=1 

Note that if a nontrivial solution to the homogeneous equations ( |2.5| ) exists, then 
by suitable scaling we can always satisfy eq. Q2.6p. We see that, at one loop, for 
d = 4 any iGraph of order 6 or higher can be decomposed in this fairly trivial way. 
A pentagon in 4-dimensions thus cannot be decomposed that way. For general d, 
iGraphs of order d + 2 or higher are decomposable. 

2.2 Reduction with coefficients linear in the loop momentum 

For a one-loop iGraph of order 5 (or lower) no trivial decomposition exists in d = 4, 
assuming that four out of the five external momenta, pj, in the loop can be considered 
as linearly independent. One way to see this is by shifting the loop momentum so that 
Vj* = 0- Then the only solution to the conditions in eq. ( |2.5|) is Xj — 0, j — 1 . . . , 5 
which is unacceptable in eq. ([2.61). We therefore turn to the next simplest possibility 
for the T's, with a linear q dependence : 

4 

Tjiq) = x j + J2 x iA ( l- t k) ■ (2-7) 
fc=i 

The single Xj is now replaced by 5 (or d + 1) variables to be determined in each 
T. Here, the four (or d) vectors must be linearly independent but are otherwise 
arbitrary. In analogy to eq. ( |2.5| ) and eq. ( |2.6| ), we now have more tensor structures 
in terms of the loop momentum : we can denote them by the shorthand 

1 , , q fl q u , q 2 q fl . (2.8) 

There are, for d = 4, therefore 1+4+10+4 = 19 independent tensor structures. Note 
that the q 2 appearing in eq. ( |2.5| ) is no longer independent since it appears as the 
trace part of q^q v . This can be extended to the inclusion of higher-rank tensors and 
other dimensions : in d dimensions, and with the inclusion of tensor up to rank k, 
we find for the number N(d, k) of independent tensor structures 

N(d ' k) = { d ~l +k ) + %{ d ~l +P ) ■ (2J) 
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In the table below we give the results for various ranks and dimensionalities. 



k 





1 


2 


3 


4 


d=l 


3 


4 


5 


6 


7 


2 


4 


8 


13 


19 


26 


3 


5 


13 


26 


45 


71 


4 


6 


19 


45 


90 


161 


5 


7 


26 


71 


161 


322 


6 


8 


34 


105 


266 


588 



Values of N(d, k) 
The number of coefficients x to be solved for is given by 



X(n, d, k) 



n 



d + k 
k 



(2.10) 



Since for d — 4 and k = 1 we have N(4, 1) = 19, it would seem that iGraphs of 
order 5 and 4 are decomposable with linear terms. However, the situation is not so 
simple since it is not obvious that the 25 coefficients for n = 5 and the 20 coefficients 
for n = 4 allow us to actually build up the 19 required tensor structures. We now 
describe how we can ascertain the number of independent structures numerically, by 
an approach that may be dubbed cancellation probing. 

We start by generating random values for the external momenta p^ and mj 
(j = 1, . . . , n). This avoids any possibility of us choosing, coincidentally, any special 
phase space point where degeneracy might occur. Then, we choose random values 
for g M precisely £ = X(n, d, k) times, and insert all this in eq. (|2.2|) . We are left with 
a set of £ linear equations for the £ unknowns x : 



(2.11) 



The £ x £ matrix M is purely numerical. We obtain it using the computer-algebra 
package MAPLE 2 which, although not numerically the fastest available, has the essen- 
tial advantage that one can easily set the precision with which numerical operations 
are performed 3 . Now, if the number of independent tensor structures that can be 
formed with our T's is less than £, the determinant of M will vanish. In an ideal 
real-number model of computation, we would thus find det(M) = 0, but in our actual 
numerical computation there will be rounding errors. A cancellation of numbers to 
'zero' will , in MAPLE, actually give a number of order 10~ p , where p is the number 
of digits specifies in the precision we tell MAPLE to use. If the matrix' determinant 



2 http:/ /www. maplesoft.com/ 
3 The relevant variable is Digits. 
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is computed by Gaussian elimination 4 , then a matrix with q zero eigenvalues will 
have a determinant of order 10~ pi3 . By letting p run down from 150 to 20 in steps of 
10, we can obtain 5 a very accurate estimate of q, especially since q must be integer. 
We give two examples to demonstrate the use of this method for the calculation of 
the zero eigenvalues. In both examples we consider a decomposition of pentagon to 
boxes, with linear and quartic terms respectively. 



The number of zero eigenvalues 
using rounding errors for the case 
of a pentagon decomposition with 
coefficients linear in the loop mo- 
mentum. The obtained value is 
q = 5.988 ± 0.098; the (abso- 
lute values of the 25 x 25 deter- 
minants range from 1.6 1 -826 to 
2.3 10- 108 . 

20 40 60 80 100 120 140 



The number of zero eigenvalues 
using rounding errors for the case 
of a pentagon decomposition with 
up to quartic in the loop mo- 
mentum coefficients. Here, q = 
188.98 ± 0.15, with the 350 x 
350 determinants ranging from 
2.7 lO" 26000 to 2.3 1 - 3320 . 



The difference £— q then gives the rank of M, and this determines the decomposability 
of the iGraph : the rank must be equal to N(d, k) for it to be decomposable. In the 
table below we give the results of cancellation probing for various n and d. 

4 This is almost unavoidable since the matrix M is not sparse, and anyway we can choose 
Gaussian elimination as an option in any case. 

5 Surprisingly, the cancellation probing appears to fail for p — 15 and p — 10, possibly since 
MAPLE may have special ways to treat these accuracies (p = 10 is default in MAPLE). 
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n 


d= 6 


d = 5 


d = 4 


d = 3 


d = 2 


d = 1 


2 


14-0 


12-0 


10-0 


8-0 


6-0 


4-0 


3 


21-1 


18-1 


15-1 


12-1 


9-1 


6-2 


4 


28 -3 


24-3 


20-3 


16-3 


12-4 


8-4 


5 


35-6 


30-6 


25-6 


20-7 


15-7 


10-6 


6 


42-10 


36-10 


30-11 


24-11 


18-10 


12-8 


7 


49-15 


42-16 


35-16 


28-15 


21-13 


14-10 


8 


56-22 


48-22 


40-21 


32-19 


24-16 


16-12 



The rank of M for various d and n, 
given as the difference £ — 

We have denoted the limit of decomposability with horizontal lines. We conclude 
that in four dimensions, n = 5 is precisely decomposable, but n = 4 is not. We now 
also see the deeper reason for this: in spite of there being 20 coefficients (one more 
than the minimum of 19), only 17 independent combinations can actually be formed. 
We also see that for sufficiently large n the number of independent combinations of 
coefficients saturates at N(d, 1) as it ought to. We conclude that in d dimensions, an 
iGraph of order d + 1 is precisely decomposable with linear terms, but one of order 
d is of course not. 

In the OPP method [18], the linear terms are precisely the spurious terms. We 
have to note that our general linear terms are not exactly those. The spurious terms 
have a specific property leading to fewer tensor structures, and we give an example 
in the appendix. Rewriting our general linear terms in terms of propagators and 
spurious terms, we see that we decompose a pentagon into boxes and triangles (like 
in [6] and [18] for example). It can be checked that the triangles always cancel, and 
therefore the decomposition is actually unique. 

At this point it must be pointed out that, in all cases where a decomposition 
is possible in principle, we actually have obtained a solution for, the system (|2.11|) . 
Once a would-be solution is found, it can easily be tested by evaluating eq. (|2.2j ) 
for additional random values of the loop momentum 6 This 'global 1=1 test' then 
verifies this solution, where of course the equality is supposed to hold only up to the 
precision used. 

This finishes the discussion for scalar one-loop iGraphs, that have unity for their 
numerator. Let us regard and iGraph of order n with a nontrivial numerator, for 
instance 

q ■ k 



D X D 2 ---D 



N 



6 In a certain sense, eq. (2.2) implies an infinite number of linear equations, of which we take £ 
and solve them. 
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Now, we can always arrange for p^ = by a suitable shift of the loop momentum, 
and write the vector as 

n 

F = ^ + 2^CiP/ , (2-12) 

3=2 

where the £'s are fixed numbers and u ■ Pj = for all j (if n > d + 1 then simply 
vanishes). We can then write 

n 

(q-k) = (q-u) + J2 Q (Dj -D 1 - fij + n x ) (2.13) 

3=2 

so that this nonscalar iGraph decomposes into scalar iGraphs of order n and 
n — 1, plus possibly a spurious term about which we do not worry since it integrates 
to zero. Our treatment of the scalar case is therefore sufficiently general. 



3. Reduction of two-loop integrands 
3.1 Preliminaries 

We now turn to the problem of reducibility at two loops. Recently several attempts 
in this direction have appeared in the literature (see i.e. [44], [40], [41]). Let us 
assume that l\ and l 2 are the two loop momenta. We consider three different kind of 
propagators for the three different loop lines of a generic two loop iGraph. 

Dih+pt), D(l 2 + Pj ), D(h + l 2 + p k ) (3.1) 

where for instance D(li + Pj) = {h + Pj) 2 — Tn 2 - and the pj are the external momenta 
associated with the propagators of the diagram. 

Such iGraphs can be denoted by the triplet (n\, n 2 , n 3 ) which indicates the num- 
ber rix of propagators that contain only the one loop momentum li, the number 
ri3 of propagators containing only the other loop momentum l 2 , and the n 2 prop- 
agators containing both. Obviously due to the symmetries of the iGraphs, for in- 
stance exchange l\ <H- l 2 , we have relations of the form (ni,n 2 ,n 3 ) <H- (n s ,n 2 ,ni) or 
(jii, n 2 , ris) -H- (jii, 713, n 2 ) provided we also exchange properly the external momenta. 
Predictably, we write the total order of the iGraph as n = n\ + n 2 + n 3 . 




h+Pa h+Pi 
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An example of the two-loop iGraph (4,2,5). 
One can see the three different loop lines. 

The propagators depending on both loop momenta are called mixed propagators. If 
these are absent the two integrals factor out and the problem becomes a double copy 
of one loop integrals. The same happens in case any other loop line is missing since, 
by shifting, one can always arrange the loop momenta such that they factor out. We 
consider these cases solved (by the one loop techniques) and will not discuss them 
further. 

The general strategy consists in finding function Xj = Xj(li,l 2 ) satisfying the 
following equation 

m 711+712 n 

J2x j D(l 1 + Pj )+ x j D(h + l 2 +p j )+ xjD(l2+Pj) = l ■ (3.2) 

j=l j=n\+l j=ni+n 2 +l 

As in the one-loop case, a generic graph of order n < 2d with 711,2,3 < 4 cannot 
be decomposed in this way, since there are Z 1)2 momenta for which all propagators 
appearing in the above equation can be simultaneously on-shell. This does not imply 
on the other hand that iGraphs of higher order must always be decomposable for any 
phase-space and mass configuration. A counterexample is the Feynman diagram of 
order 5 in two dimensions 7 : if all internal lines in this self-energy Feynman diagram 
are massless, it is possible to choose the two loop momenta components such that 
all five propagators are simultaneously on-shell. 




The requirement for trivial decomposition (with coefficients Xj that are constants 
with respect to the loop momenta) now reads 

ni ni+n-2 n 

J2x j D(l 1 +p j )+ x j D(l 1 + l 2 +p j )+ XjD(l 2 + Pj ) = l . (3.3) 

j = l j = m + l .3=711+712 + 1 

Proceeding in analogy with our one-loop discussion, we find that we have to satisfy 
the following equations : 

ni+n 2 n rti+n 2 

Y x i = Yl x i = Yl x i = ' ( 3 - 4 ) 
j=i j=m+i j= ni +i 

7 This diagram can, in fact, be decomposed, but not by the method described above: instead one 
has to use IBP techniques. 
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ni+n 2 n 
j=l i=ni+l 

and 

n 

In total there are 2d + 4 conditions, so that the minimum size of a trivially de- 
composable iGraph is 2d + 4. In four dimensions, scalar iGraphs can therefore be 
decomposed down to n = 11. Again in analogy, for n = 11, since by shifting we can 
arrange pi + • • • + p ni +n 2 — as well as p ni +i + ■ ■ ■ + Pn = 0, the only solution to 
the 2d + 3 homogeneous equations is Xj — 0, j — 1, . . . , n, and this fails the inho- 
mogeneous equation. On the other hand, since any subset of an iGraph is itself an 
iGraph, any iGraph with ri\ > 6, n 2 > 6, or n 3 > 6 is trivially decomposable (for 
d — 4). Furthermore, with linear terms we see, from the one-loop discussion 8 that 
we only have to consider two-loop iGraphs with 

"1,2,3 < 4 (= d) , n 1 + n 2 + n 3 <ll(=2d + 3). (3.7) 

A word of caution is in order here. We may have a case where n\ + n 2 > 6 or 
"2 + "3 > 6 and then decide to perform a decomposition a la one-loop with trivial 
terms, taking Zi as the loop momentum, for instance, and l 2 as one of the "external" 
momenta or vice-versa. Since the solution of the eq. ( |3.3| ) is nonlinear with respect 
to the external momenta, due to matrix inversion, the resultant decomposition will 
not have, in general, the simple form of iGraphs again, and the emerging integrals 
will belong to very different classes of functions. 

The number of two-loop iGraphs that we have to consider is therefore not very 
large : 4 for d = 2, 10 for d = 3, 19 for d = 4. 

3.2 Further reduction with linear terms 

With trivial decomposition we see that we can always end up with an iGraph of 
order 2d + 3. Like in the one loop case, we now add coefficients linear in the loop 
momentum and hope for further reductions. 

A note is in order here. In the one loop case the resulting integrals were always 
scalar. The reason is that any contraction of the loop momentum with any vector 
can either reconstruct denominators or be a spurious term. After integrating, in the 

8 In case there at least 6 propagators in one loop line we can first reduce the propagators in this 
loop line with constant coefficients and then continue further if possible. In the case of 5 propagators 
we already know that adding linear terms that depend only on the loop momentum of this loop 
line and take all coefficients that depend on the other loop momentum to zero, we can again solve 
the problem a la one loop. 
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case a denominator is reconstructed the remaining integral is a scalar integral with 
fewer denominators. In case the term is spurious it vanishes after integration. In 
two loops this is not the case anymore. One can always use dot products of the loop 
momenta with the momenta of the integrals to write relations like 



The denominator D(li + pj) may, however, not be present in the integral in case 
Pj appears in a propagator of another type such as D(l 2 +Pj) 2 . Then, the product 
li ■ pj may be an irreducible scalar product (ISP) [40]. But not always; the ISP's 
of an integral are more complicated to write. For example, if there are enough 
propagators of the type (Zi +pi) 2 in the diagram such that the p;'s can form a basis, 
one can rewrite linear combination of the p^s and manage to reconstruct 

denominators. There is a specific number of ISP's in any diagram and one can have 
some freedom in how to write them. The integrals of the resulting basis can have 
numerators with ISP's in some power. It is not obvious that a scalar integral with 
a specific number of denominators is more difficult to calculate than a non-scalar 
integral with fewer denominators; however it is commonly accepted that this is the 
case. 

Again, we want to write, if possible, the number 1 as in eq. ( |3.2|) . We use general 
linear terms in the sense that in every dimension we construct a basis ti (possibly, 
but not necessarily, the external momenta in the iGraph) and we have 



with the ai,bij, and c^- constants with respect to the loop momenta. Since in d 
dimensions, we need d vectors to construct such a basis, it is obvious that for an 
iGraph of order n we start with (2d + l)n coefficients. As in the one-loop case, we 
give a table that contains, for every dimension we worked with, the number of tensor 
structures (denoted by T(d)) and the number of independent coefficients that we 
have explicitly calculated by cancellation probing as described above. 



2h ■ p 3 = [(h + Pj ) 2 



2l ;2 2,2 

rrijl — l 1 — pj + rrij 





(3.9) 



T(d) = (4d 2 + 18d + 2)/2 
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n 


d = 6 


d = 5 


d = A 


d = 3 


d = 2 


d = 1 


3 


39-0 


33-0 


27-0 


21-0 


15-0 


9-0 


4 


52-0 


44-0 


36-0 


28-0 


20-0 


12-2 


5 


65-1 


55-1 


45-1 


35-1 


25-1 


15-5 


6 


78-3 


66-3 


54-3 


42-3 


30-3 




7 


91-6 


77-6 


63-6 


49-6 


35-8 




8 


104-10 


88-10 


72-10 


56-10 






9 


111-15 


99-15 


81-15 


63-17 






10 


130-21 


110-21 


90-21 








11 


143-28 


121-28 


99-30 








12 


156-36 


132-36 










13 


169-45 


143-47 










14 


182-55 












15 


195-55 












T(d) 


127 


96 


69 


46 


27 


10 



In the table a line distinguishes between reducible and non-reducible cases. Re- 
ducibility is explicitly checked by eq. (|3.2p. For all iGraphs defined by eq. (|3.7|), 
we see that the number of independent coefficients becomes equal to the number of 
independent tensor structures, when reducibility is attained. In all dimensions every 
n = 2d + 2 case is reducible with linear terms to a n = 2d + 1 iGraph. In four 
dimensions, we can decompose every integral down to integrals with 9 denominators. 
At this point, we are one step away from the limiting value of 8. 

3.3 Comments on the number of independent coefficients 

The most difficult part in our counting is always the number of independent co- 
efficients. As shown, we do it numerically but we would like to understand more 
the reason why we have as many independent coefficients as we do find. We try to 
demonstrate here a way to estimate this number; for the case of linear terms we will 
give some examples. 

We can rewrite the terms in eq. ( |3.9| ). As we mentioned, terms of the form 
li ■ pj either reconstruct denominators either become ISP. Let us assume the (4, 1,4) 
iGraph in four dimensions. It has in principle two ISP that we call <j\ and Oi. To see 
this, we note that in any of the loop lines that consists of four denominators, there 
are three external momenta. We can always "borrow" a fourth one from the other 
line to have a complete basis and write any product k linear combination 

of the four propagators and an ISP. The ISP in that case would be the product of U 
with the momentum we borrowed. Repeating for the other loop line we get the two 
ISP's. Using the coefficients of eq. ( |3.9| ) in eq. ( |3.2| ) we can either get a denominator 
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times a constant or an ISP, or products of two denominators. We write this equation 
schematically as 

9 9 9 

i = e a{i, + E E (3.10) 

i=l i=l j,j~>i 

This means that by writing {1, a%, o^Di there is a constant coefficient in front of 
every of these different terms: 

{1, 0i, cr 2 } = a • 1 + &i<7i + & 2 0"2 

where a, 61, 62 general numbers. In our particular example, we have 9 x 3+45 x 1 = 72 
coefficients. However, we started with a problem with up to cubic power in loop 
momenta and ended with up to quartic powers since we have these products of two 
denominators times some constants. These higher powers have to cancel, which 
means that we have to put extra constraints on our coefficients. We have 6 such 
constraints to cancel, namely the 

^2? 4' (A ' (h ' h)^, (h ' 

terms. As a result we end up with 72 — 6 = 66 independent coefficients, which is the 
number we get numerically as well. If we now try to decompose an iGraph of order 
10 we can prove that eq. ( ETTDj ) becomes 

10 99 

1 = e a(i, + E E D * D >w ( 3 - n ) 

i=l i=l j,j>i 

We don't need to go up to 10 in the product of 2 denominators since 9 

D 10 oc (Di, ...,D 9 ,ai,a 2 ) 

We still need them, though, in the first term to produce terms of the type (<Jj) 2 . In 
that case we have 69 independent coefficients and this graph is reducible. Adding 
more propagators we do not get more independent coefficients. In the same way 
one can count the independent coefficients in all dimensions although it is clear that 
it is safer to find their number numerically since there are a lot of overlaps in the 
tensor structures for higher cases. The way of rewriting the general linear terms 
as propagators and ISP's in the example above is still not the OPP method. In 
an extension of the OPP method to two loops, one would find the ISP's of every 
subdiagram and would avoid terms like Df. We expect something similar to the 
one-loop case to happen then, rewriting o\ and 02 in the form of true ISP's of every 
subdiagram's contributions of the terms with the highest number of denominators to 
cancel. For that it is possible that special properties exist, as again in the one-loop 
case where spurious term solve a lot of equations by putting automatically tensor 
structures to zero. 

9 This is actually the point where the number of dimensions plays a role in the counting. 
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3.4 Reduction with general quadratic terms 

We try to reduce our iGraphs further with the use of general quadratic terms in the 
coefficients. In this case the coefficients become 

Xj = eg + } J bij{li ■ tj) + ^2 °y ■ tj) + ^2 dijkih ■ tj)(h • h) 

j j j<k 

+ e ijk(h ' tj)(h ■ tk) + fijh(h " tj)(h ■ tk)) 

3<k j,k 

(3.12) 

We give the number of tensor structures T(d) and the original number of coefficients 
with general quadratic terms, Ci(d). The coefficients depend on the number of 
propagators n. Notice that the expression for T(d) is not valid for d = 2. In that 
case, there is more overlap between the highest tensor structures. More specifically, 
for this particular case one can completely reconstruct the Ifl^l^ structure from l\ 
and (/i • lz)l]Vi allowing fewer independent structures to be constructed. 

T(d) = 4d 3 /3 + I0d 2 + 20d/3 - 2 (3.13) 
C^d) = (2d 2 + 3d + l)n (3.14) 

With quadratic terms we start with (2d 2 + 3d + 1) x n coefficients in total, not all 
of them of course being independent. Using cancellation probing we are able to find 
the number of independent coefficients for different iGraphs in different dimensions. 
We put our findings in the following table : 



n 


d = 4 


d = 3 


d = 2 


3 


135-4 


84-3 


45-3 


4 


180-6 


128-6 


60-6 


5 


225-18 


140-16 


75-15 


6 


270-38 


168-32 


90-30 


7 


315-65 


196-53 




8 


360-98 


224-80 




9 


405-136 


252-108 




10 


450-180 






11 


495-225 






T(d) 


270 


144 


60 



In the table above a line distinguishes again between reducible and non-reducible 
cases, according to eq. ( |3.2|) . As we see, we can decompose in 2 dimensions an iGraph 
of order 5 to lower order iGraphs as indicated by unitarity. However, there is no such 
solution in 3 and 4 dimensions and in this case we have to investigate what happens 
if we add cubic, quartic terms and so on. The two-dimensional case is exceptional 
here because of the extra properties that lower the number of tensor structures. 
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3.5 Reduction with general cubic terms 

We focus now on 3 and 4 dimensions since we finished the reduction in 2 dimensions. 
We include general cubic terms and our coefficients become 

Xi = di +y^ j b ij (i 1 ■ tj) + y^cj 3 -(z 2 • tj) +^2fd ijk (ii ■ ■ t k ) 

j j i<k 

+ ^ Zijk(h ■ tj){h ■ tk) + ^ fijk(h ' tj){h ' tk) 
j<k j,k 

+ ^2 9ijki(h ■ tj)(h ■ *fc)(Zi • u) + /lijfcK^i ■ ■ *fc)(^ • *o 
j<fc<« j j<fc 

2 i<fc j<fe<2 

(3.15) 

We give the number of tensor structures T(d) and the original number of coeffi- 
cients with general cubic terms C\{d). The coefficients depend on the number of 
propagators n. The expression for T(d) is valid for d > 3. 

T(d) = 2d A /3 + 22d 3 /3 + 71rf 2 /6 + d/6 + 1 (3. 16) 

d(d) = (4d 3 /3 + 4d 2 + lld/3 + l)n (3.17) 

This means that in d dimensions we start with C±(d) coefficients, not all of them 
being independent. We run the MAPLE code in the case of the iGraph of order 7 in 
3 dimensions and we find that out of 588 original coefficients, 360 are independent. 
This is the number of tensor structures as well. Using another PYTHON-based 
program 10 , we can actually solve the system decomposing any iGraph of order 7 in 
3 dimensions to lower iGraphs with general cubic terms, and perform the 1=1 test. 
This means that with cubic terms we are able to decompose any two-loop iGraph 
in 3 dimensions to up to a 2d iGraph as expected from unitarity. In the same way, 
we can investigate d — 4, and we get a valid decomposition: from our original 1485 
coefficients, 831 are independent and all the tensor structures can be reconstructed. 
Actually we did the same in 5, 6, 7, and 8 dimensions, and managed to decompose 
every integral of generic order 2d + 1 or higher, using cubic terms. We believe that 
this is a general result for any dimension, except of course for d = 2, where the same 
can be achieved with only quadratic terms. 

4. Conclusions 

By introducing the notion of iGraphs, we have investigated the decomposability of 
Feynman amplitudes at one and two loops. In both cases we have demonstrated (in 

10 http://codepad.org/zT4wUxCJ 
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the two-loop case, up to d = 8) that generic iGraphs can be decomposed down to 
the unitarity-based limit of n — d and n = 2d, respectively. For one-loop graphs, the 
inclusion of linear terms is sufficient, and we have elucidated the relation between 
these linear terms and the OPP's spurious terms. At the two-loop level, ultimately 
cubic terms are needed (quadratic for d — 2), and the decomposition is seen to 
lead to non-scalar, non-vanishing integrals. If one wants to design a two-loop OPP 
method, it is clear that one has to take our general linear, quadratic, cubic terms 
and rewrite them in terms of propagators (this would lead to contributions with less 
denominators) and ISP's of each subdiagram seperately 11 . Our work is basically the 
starting point of an OPP method and a proof that reductions that were conjectured 
in [40, 41] are actually valid and survive the global 1 = 1 test, in case one includes 
all relevant cuts. 

We note once again that the resulting integral basis, obtained in this way, is 
clearly not a minimal one. We are aware of cases that could be further decomposed 
using IBP identities. Achieving a proper level of understanding of the interplay be- 
tween OPP and IBP will certainly open the road for an efficient reduction of two-loop 
amplitudes at the integrand level. 

Acknowledgement: We acknowledge useful discussions with Prof. M. Czakon. 

A. Appendix: Pentagons to Boxes with the OPP method 

By defining 

D n ( Pl ,p 2 , ...,p n ) = — — — 1 (A.l) 

D(q + p 1 )D(q + p 2 )-D(q + p n ) 

we try to decompose this pentagon 

to the 5 following boxes 

Da(P2,Ps, Pa, Pb), D 4 (P3, Pa, Pb, Pi), D a (pa, Pb, P\,Pi) 

Da(P5,Pi,P2,P3)Da(pi,P2,P3,Pa) 

By shifting the loop momenta we can always arrange that the momenta pi of the 
denominators sum up to zero: 



i=i 



ii 



In the same way that one-loop OPP has different spurious terms for every subdiagram 
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We write the spurious terms in the following form 



Si = e(q + p2,q + P3,q + P4,q + p 5 ) = q^(p 3 , p 4 , p 5 ) - e M (p 2 ,P4,P5) + 
^(P2,P3,Ps) - ^(p 2 ,p 3 ,p 4 )) + e(p 2l p 3} p 4 ,p 5 ) = 
+ B x 

(A.2) 

We then have: 



K = ^(j>3,Pi,P$) ~ ^(j>2,Pi,P$) + ^(P2,P3,P5) - ^(P2,P3,Pi) 

K = e M (P4,P5,Pi) - ^(P3,P5,Pl) + e M (P3,P4,Pi) - e M (p 3 ,P4,P5) 
A 3 = ^(P5,Pi,P2) - ^(p4,pi,p 2 ) + e^(p 4 ,p 5 ,p 2 ) - e^(p 4 ,p 5 ,Pi) 
K = e^(p l ,p 2 ,p 3 ) - e tM (p 5 ,p 2 ,p 3 ) + e"(p5,Pi,p 3 ) - e"(P5,Pi,P2) 

^5 = e>1 (P2,P3,P4) - ^(Pl,P3,PA) + ^(PuP2,Pi) ~ ^{p U p 2 ,Pz) 

Using the OPP master formula [18] we get 



(A.3) 



1 = [tZi + diSi]£)(g + pi) + [d 2 + d 2 5 2 ] J D(g +p 2 ) + ... + 
[d 5 + d 5 S , 5 ]D(g + p5) 



(A.4) 



We compare the polynomials of the two sides of eq. ( |A.4 ). The q 2 q ll terms must 
vanish for any value of q which means 



£(*4*) = 



(A.5) 



Notice that from eq. ( A.3 ) 



i>=o 



»=i 



which means that if all di are equal eq. (|A.5|) is satisfied. 

This is actually the only solution. To see that, consider the sum 



<M£ + • • • + d~ 5 ^ 



Substituting 
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Pb = -Pi -P2-P3-P4 



we get 



J2(di A i) = (*i + ^2 - 4d 3 + di + d 5 Vij>i,P2,Pi) 
+ (-Ji + 4d 2 - da - ^4 - d 5 )e M (pi,p 3 ,p 4 ) 
+ (-di - d 2 - d 3 + 4J4 - d 5 )e M (pi,p 2 ,p 3 ) 
+ (-4di + d 2 + d 3 + d 4 + d5)^(p2,P3,P^) 

(A.6) 

In 4 dimensions, for this sum to be zero, all coefficients in front of the vectors 
are zero (the vectors are linearly independent). This results in the system: 



d\ + d 2 — 4d 3 + d 4 + d 5 = 

—d\ + 4d 2 — d 3 — d 4 — G?5 = 

— c?i — d 2 — d 3 + 4d 4 — d 5 = 
—4c?! + d 2 + d 3 + d 4 + ds = 

(A.7) 

One can solve the system and find 

di = d 2 = d 3 = d 4 = ds (A. 8) 

The only solution is when all coefficients are equal and from now on we call them 
d. Then we take a look at the quadratic in q parts of the right-hand side of eq. ( A.4j) . 



They cancel as well but one has to be careful since they come from two terms, the 
q 2 and the q^q v term. We look at the latter term 



dq^^{q -Pi) A? 



i=l 



We have : 



i=l 



(A 2 - A 5 )(q ■ p 2 ) + (A 3 - A 5 )(q ■ p 3 ) + (A 4 - A 5 )(q ■ p 4 )) 

(A.9) 
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where p§ = 



—px — P2 — P3 — P4 is again used. With this substitution we have: 



Ax- A 5 = -5e' 1 {p 2 ,p 3 ,p 4 ) 
A 2 -A 5 = 5e"(pi,p 3 ,p 4 ) 
A 3 - A 5 = -5e' 1 (px,P2,PA) 
A 4 -A 5 = 5e^{p 1 ,p 2 ,p 3 ) 

(A.10) 



Inserting the Schouten Identity 



e(Pl,P2,P3,P4)<7 M = ^{P2,P3,P4){q ■ Pi) - ^{Pl,P3,P4,){q ■ P2) + 
^(Pl,P2,P4)(q-P3) ~ £ fl (Pl,P2,P3)(q-P4) 



(A.ll) 



we get 



d q^^2(q-Pi) A f = -5dq 2 e(px,P2,P3,P4) ( A - 12 ) 

i=l 

That is exactly the property that the spurious terms have that makes the solution 
to the system possible. The q^q v terms all vanish owing to the Schouten identity 
except for the trace part proportional to q 2 solving 9 out of 10 equations in one go. 
The total number of nontrivial equations is therefore 10 and not 19 in this case and 
the system has a solution. To complete the story, taking into account eq. ( [A. 8 



we 



are left with 5 di and one d, which are uniquely now determined from the remaining 
6 equations, namely constant part, q 2 term and q^ term. 
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